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Abstract 

We use the recent estimates of NNLO splitting functions, made by van Neerven and 
Vogt, to perform exploratory fits to deep inelastic and related hard scattering data. We 
investigate the hierarchy of parton distributions obtained at LO, NLO and NNLO, and, 
more important, the stability of the resulting predictions for physical observables. We 
use the longitudinal structure function Fl and the cross sections aw , oz for W and Z 
hadroproduction as examples. For Fl we find relatively poor convergence, with increasing 
order, at small x; whereas crw,z are much more reliably predicted. 



1 Introduction 

With the increased precision of deep inelastic scattering data JTJ , and the need for accurate pre- 
dictions at the Tevatron and the LHC, it is clearly essential to extend global parton analyses 
to next-to-next-leading-order (NNLO) in as- Although the relevant deep inelastic coefficient 
functions have been known for some time [[|, there is only partial information on the corre- 
sponding splitting functions. The N = 2,4,6,8 (and 10 for non-singlet) moments have been 
calculated ||, which effectively provide information on the high x behaviour of the splitting 
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functions. Also known is the most singular logl/x behaviour at small x, both for the singlet 
[|J, and the phenomenologically less important nonsinglet || splitting functions, and the lead- 
ing nf contributions |J of the nonsinglet splitting functions, and of the Ca dependent part of 
P gg [0. Recently van Neerven and Vogt M have constructed compact analytic expressions for 
the splitting functions which represent the fastest and the slowest evolution that is consistent 
with the above information. We believe that these two extreme behaviours are indeed realistic. 
Although there are indications that the true behaviour of the splitting functions is likely to be 
slightly nearer to that corresponding to the slow evolution possibility^, for simplicity we shall 
use the average of the two extremes for our 'central' NNLO analysis. 

It is important to stress an important difference between our analysis and the procedure 
used by van Neerven and Vogt || . The latter authors start from a fixed set of partons and 
a fixed scale (~ 30 GeV 2 i.e as = 0.2) and present the differences between LO, NLO and 
NNLO evolution. Here we compare the partons, and the consequent predictions for physical 
observables, obtained by performing global analyses at LO, NLO and NNLO. Both works 
present NNLO results obtained using the extreme estimates of the 0(a s ) splitting functions. 

In Section 2 we discuss the changes to the global analysis that are necessary in going from a 
NLO to NNLO formulation. Then, in Section 3, we present seven new fits to the deep inelastic 
and related data; that is LO, NLO and five NNLO analyses. To gain insight into the impact 
of the NNLO contributions, we discuss essential features of the fits in terms of the behaviour 
of the splitting (and coefficient) functions. In Section 4 we compare the partons obtained in 
the LO, NLO and NNLO analyses, paying particular attention to the gluon distribution in 
the small x region. The parton distributions are scheme dependent and are not themselves 
observable. The comparison of LO, NLO and NNLO predictions for physical observables is 
much more meaningful. In Section 5 we study the predictions for the longitudinal structure 
function, F^. This is a particularly relevant observable as it directly reflects the behaviour of 
the gluon distribution at small x, and hence most directly probes the stability, or convergence, 
of parton analyses as we go from the LO, to the NLO, and then to the NNLO framework. In 
Section 6 we compare the LO, NLO, NNLO predictions for the cross sections of W and Z boson 
production at the Tevatron pp collider and at the LHC These observables mainly depend on 
the quark distributions in the region Q 2 ~ 10 4 GeV 2 , and x ~ 0.05 and 0.006 respectively. The 
stability of the predictions offers the possibility of using the W and Z events as a luminosity 
monitor of the collider. Finally in Section 7 we give our conclusions. 



2 Global analyses at NNLO 

The procedure is based on the NLO analyses described in Refs. [§, |i~0"| . However at NNLO it 
is important to allow the gluon distribution to become negative in the low x, low Q 2 domain. 

1 A view confirmed by private communication with A. Vogt. 
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We therefore adopt the parameterization 

xg{x, Ql) = A g x~ x ° (1 - xfo (1 + e a Vx~ + l a x) - A' g x~ x * (1 - xf'» (1) 

at the starting scale Q 2 = Ql = 1 GeV 2 of the evolution. The parameter rj' turns out to be 
large in the additional negative term and so this contribution is only important at small x. 

It is necessary to implement other extensions of the formalism when going to NNLO. First, 
we use the three-loop expression for as, in the MS scheme. Second, we require more detailed 
matching conditions when evolving through the heavy flavour thresholds. The NNLO treatment 
of heavy flavours is discussed in the Appendix. 

Our main interest is in the quality of the fit to deep inelastic data at small x. At high x 
we have a slight inconsistency in our NNLO analyses in that we use NLO expressions to fit to 
Drell-Yan, jet production and boson rapidity asymmetry. The NNLO corrections to all 
these quantities have not yet been calculated. However note that the physical observables that 
we study (namely F L and (Jw,z) sample low x partons, which are determined mainly by deep 
inelastic data for which the NNLO formalism is consistent. 



3 The new global fits 

We perform LO, NLO and NNLO global fits to the set of deep inelastic and related data that was 



used in Refs. [|| fLOl , except that now we use the jet Et distribution measured at the Tevatron 
to pin down the gluon distribution at large x, instead of prompt photon hadroproduction. The 
QCD description of the latter process has outstanding theoretical problems JTTJ]. A second 
change is that we include all the available preliminary HERA data Pj], which have higher 
precision than hitherto. 

The consequence of replacing prompt photon data by the jet data is that the NLO fit is 
now similar to that achieved by the previous MRST(g|) set of partons |], |Kj. A satisfactory 
description of the Tevatron jet data is obtained, including particularly the normalization. 

Five NNLO fits were performed. The 'central' fit and the four extremes (A q A g , A q B g , B q A g , 
B q B g ), where A^Bi) corresponds to the slow (fast) evolution of parton i — q,g. It turns out 
that the NNLO fits with slow and fast gluon evolution are very similar, and so it is sufficient 
to present results for just two of the extreme choices of the splitting functions, namely 

A = A q Ag (slow evolution), 

B = B q B g (fast evolution). (2) 
In Figs. 1 and 2 we show the LO, NLO and NNLO descriptions of the F 2 data [12| in a 



few representative x bins. We display only the 'central' NNLO fit. However the quality of all 
the NNLO fits is similar. It is encouraging to note that, as we proceed from the LO — > NLO 



3 



— > NNLO analysis, there is sequential improvement in the overall quality of the description of 
the data. In particular, in going from the NLO — > NNLO fit, there is an improvement in the 
simultaneous description of the NMC and HERA F 2 data. Indeed the quality of the NNLO fit 
is improved for almost all subsets of the data. 

From Figs. 1 and 2 we can see that at NNLO the scaling violations increase both at small 
x and at large x. At small x this is due mainly to the NNLO contribution to P qg , whereas at 
large x the NNLO term in the coefficient function plays the dominant role. The relevant x —>■ 
behaviour of the splitting functions are| 

. . as 1 / 5«5 14a%hi(l/x) \ , . 

P aa (x) = 2n f —-[l + - + — K - J -^- + ---) 3 

qgK ' / 27r3V 3a; 9 x ) v; 



Pgg \ x ) 



X 



1 - 



61n/a s /395 1 llvr 2 / 295 _ vr 2 \ \ a|ln(l/x 

36 ~~x ~ V104 ~ 2^ ^ ~ ~72 + n/ V2808 ~ 972// x 



,(4) 



where as = (3/n)as, and the x — > 1 behaviour of the quark contribution to the F 2 coefficient 
function is 

^ / x r/, \ as „ /41n(l — x) — 3 
C 2 , q (x) = 5(l-x) + ^C F l K K — 



x 



+ 



An) 



.(5) 



As well as the improvement in the quality of the fit, we can investigate the importance of 
the increased scaling violations by looking at the higher-twist component of F 2 extracted using 
a phenomenological analysis in which a term (D(x)/Q 2 )F 2 (x,Q 2 ) is included in the fit, as in 
Ref. fL3| |. The values of the higher-twist coefficient D(x) can be seen in Table 1. At very high x 



a large positive higher-twist contribution is clearly needed. This decreases slightly as we move 
from LO to NLO to NNLO, but there is no indication that its presence will be eliminated by 
even higher orders. We note that the conclusion that NNLO contributions largely remove the 



need for higher twist at high x in previous NNLO analyses [T3] has been based on analysis of 
CCFR data only, which exists at far higher W 2 than the SLAC data included in our higher- 
twist fit, though it has also been suggested that when NNLO coefficient functions are used the 
higher twist may be almost entirely due to target mass effects [0. At x — 0.4 — > 0.5 the 
higher-twist contribution changes sign, becoming generally negative. At LO its magnitude is 
then quite large, demonstrating that the evolution is too slow at low Q 2 , both for NMC and 
HERA data, as is obvious from Fig. 1. The magnitude of the higher-twist contribution for 
x < 0.3 decreases significantly going to NLO, and decreases again, to very small values, at 
NNLO. Indeed, the sign of the small-x higher-twist contributions at NNLO is not even well- 
determined, with many x-bins preferring a slightly positive value. The implication seems to be 
that higher-twist contributions at small x are small, and their apparent size is decreased by the 
inclusion of more perturbative corrections. 

2 For the LO splitting function Pqg\x) we use the coefficient of the moment space expression in the limit 
N — > rather than the real limit as x — > 0. 
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Table 1: Values of the higher-twist coefficient D(x) extracted from the LO, NLO and NNLO 
fits 



X 


LU 


ATT ( \ 


ATATT (~\ 


- 0.0005 


f\ /I T r A 

—0.4754 


0.0116 


—0.0061 


0.0005 - 0.005 


—0.2512 


—0.0475 


0.0437 


0.005 - 0.01 


—0.2481 


— 0.1376 


n r\r\ a o 

—0.0048 


0.01 - 0.06 


—0.2306 


— 0.1271 


—0.0359 


0.06 - 0.1 


r\ 1 o to 

—0.1373 


r\ AO01 

—0.0321 


0.0167 
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— U. izoo 


— U.Uooi 


U.UU ( 


0.2 - 0.3 


-0.1210 


-0.0893 


-0.0201 


0.3 - 0.4 


-0.0909 


-0.1710 


-0.1170 


0.4 - 0.5 


0.1788 


-0.0804 


-0.0782 


0.5 - 0.6 


0.8329 


0.3056 


0.1936 


0.6 - 0.7 


2.544 


1.621 


1.263 


0.7-0.8 


6.914 


5.468 


4.557 


0.8 - 0.9 


19.92 


18.03 


15.38 



For each fit - LO, NLO, NNLO - we use the one-, two-, three-loop expression for the (3 
function, e.g. in the NNLO fits the connection between as and Aj^g involves fa evaluated in 
the MS scheme. For completeness we show in Table 2 the values of the QCD coupling, together 
with Aqcd, found in the different global fits. For the NNLO fits the value of Aqcd is kept the 
same for the extremes as for the central fit, but would change by only a tiny amount if left 
free. In the fits where a higher-twist component is allowed, at each order the extracted value 
of Aq CD increases, reflecting the effect of the increased scaling violation by the new data at low 
Q 2 included in these fits. This increase in Aqcd is only 10 — 15% (decreasing with increasing 
order), leading to a corresponding increase in a 5 (M|) of 0.002 — 0.003. 

4 Implications for parton distributions 

In Fig. 3 we compare the parton distributions found in the NNLO fit to those in the NLO anal- 
ysis. We plot the NNLO/NLO ratios for the gluon, and the up and down quark distributions, 
at two values of Q 2 . 

As we go from the NLO to the NNLO analysis, several changes in the distributions are 
worth noting. First, the decrease of the quark distribution at high x and the slight increase at 
low x reflect the behaviour of the coefficient functions Ci, q and Ci, g respectively. Second, recall 
that in the NLO analysis the input gluon distribution decreased at small x. At NNLO we see 
the gluon decreases even more. This decrease at low x occurs because of the increase of P qg , see 
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Table 2: The QCD coupling and A parameter 





a s {M z ) 


(n /= 4) /n„Trt 
A LO or MS ( GeV ) 


LO 


0.1253 


0.174 


NLO 


0.1175 


0.300 


(NNLO) central 


0.1161 


0.242 


(NNLO)a 


0.1161 


0.242 


(NNLO)b 


0.1161 


0.242 



(IH). The consequent rise at x ~ 0.1 is to ensure that the momentum sum rule is satisfied. The 
gluon distribution drives the evolution at small x. As we evolve to higher Q 2 , the effect of the 
NNLO term in the splitting function decreases, and a smaller gluon leads to slower evolution 
than in the NLO analysis. Hence, for example, by Q 2 ~ 10 4 GeV and x ~ 10 -4 , all NNLO 
partons are about 10-15% smaller than those at NLO. 

Since the biggest NNLO effect is in the small x behaviour of the gluon, we study this 
distribution in more detail. Fig. 4 shows the gluon obtained in the LO, NLO and NNLO global 
fits at various values of Q 2 . A clear LO — > NLO — > NNLO hierarchy^of the small x behaviour 
of the gluon is evident, which reflects the direct link with the HERA deep inelastic data via 
P qg of (|3|). Note also that the evolution of the NNLO gluon is made even slower because of the 
(small) negative NNLO contribution in P gg , see (f|). 

The 'starting' parametric forms of the gluon found in the LO, NLO and NNLO global 
analyses are 

5.23v/i+7.33x) (LO) 
3.92v^ + 4.68x) - 1.67x-° 032 (l-x) 8 - 21 (NLO) (6) 
3.22v/i + 4.36x) - 0.705a;- ai51 (l - x) 8 ' 69 . (NNLO) 

The 'extreme' curves A and B, plotted in Fig. 4, demonstrate that the greatest uncertainty, 
coming from the lack of complete knowledge of the NNLO splitting functions, is in the small x 
behaviour of the gluon. Nevertheless even allowing for the 'extreme' spread in the NNLO fits 
we see that the hierarchy in the small x behaviour of the gluon persists. 

Fig. 4 also shows that the gluon obtained from the NNLO analysis becomes negative at small 
x and small Q 2 , as anticipated in (|l|). However the gluon distribution itself is not a physically 
observable quantity. It is scheme dependent. For example, Fig. 4 shows the gluons obtained 

3 The wobble seen in the LO gluon at x ~ 0.1 for Q 2 < 20 GeV 2 is a consequence of momentum conservation 
and a much too large a gluon at small x. 



xg{x,Q ) = 




31.2^.390^ _ ^6.18 (!_ 

51.8x°' 535 (l - x) 6 - 55 (1 - 
14.4x a397 (l - s) 5 - 56 (1 - 
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from analyses at different orders in the MS factorization scheme. If on the other hand we were 
to adopt the DIS scheme, then we find that the NNLO gluon is only marginally negative at 
low x at Q 2 = 2 GeV 2 . In order to investigate the true implications of the convergence of the 
perturbative series we must examine the predictions for physically observable quantities. The 
behaviour of the longitudinal structure function, Fl, is particularly appropriate as it is sensitive 
to the small x behaviour of the gluon. The production cross sections of W and Z bosons at 
the hadron colliders are representative of other relevant observables. We therefore study the 
predictions for these quantities below. 



5 Predictions for Fl 



The LO contribution to Fl is 0{as), and so a consistent (factorization scheme independent) 
NNLO prediction of Fl requires the 0(a\) coefficient functions. These are not known at 
present, but we do know much of the same information as for the 0(ag) splitting functions, 
that is the N = 2, 4, 6, 8 moments and the x — > behaviour. Hence we estimate the coefficient 
functions in the same spirit as used by van Neerven and Vogt for the 0(o? s ) splitting functions. 
The 'central' estimates for the NNLO contributions to Cl are (where the common factor of 
(as/(47r)) 3 is taken out) 



C { S 9 {x) 



ln(l/x) 



1200 



x 



x 



+ 1095 ln 2 (l/x) - 5960 + 21512x 2 + 1928 ln(l 



C 



(3) 
L,PS 



» 



n f (^381 

+n 2 f (— + 148.8 ln 2 (l/s) - 5 - 74Lr 2 - 1471n(l - x)^j 

-323 ln 2 (l/x) - 3916 - 47526s 2 - 21954 ln(l 
+n f (-89 ln 2 (l/x) + 863 + 2796x 2 + 2038 ln(l - 
+n^l51n 2 (l/x) - 54.3 + 72.4x 2 - 231n(l 

n f( 



x 



(7) 



(8) 



( 169 ln(l/x) 700 + lgg x ^^ jx) + ^2 + 42 g ln(1 _ %) + 3m \ 



x 



+n) (— + 61 ln 2 (l/a;) - 25 + 42.7x 2 + 7.2 ln(l - x)j 



(9) 



where NS and PS refer to quark non-singlet and pure-singlet respectively. In fact the n 2 
dependent part of the non-singlet coefficient function is in principle known exactly from the 



calculations in |16[], but are small and well modelled by our simple analytic expression. 

The behaviour of the Fl gluon coefficient function is shown in Fig. 5. The two dominant 
features are (i) a sizeable contribution just below x — 1, and (ii) a large growth with decreasing 
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x arising from the most singular terms found in Ref. 19. In fact at small x we have[] 

c -m^Ht + (t-« 2 >)^] ™ 

and the same expression, modulo the colour factor Cf/Ca = 4/9, for Cl, ps (except at leading 
order) . 

The non-singlet coefficient functions beyond LO are very strongly peaked as x — > 1. At 
NLO the coefficient function [pj]] (with (a,g/(47r)) 2 factored out) in this limit behaves like 



2C F ln 2 (l - x) + (9 - 8C(2))CVln(l - x) 

(11) 



+ (4(C(2)-l)^-(ll-^ / ))ln(l-x 

and there is an enhancement compared to the LO result, (a,g/(47r))4CV£, due to the ln(l — x) 
terms. The machinery for computing the dominant ln(l — x) terms for Cl,hs for all orders 
in as has recently been devised fTS], and in principle we could use this to evaluate the parts 



oc ln m (l — x) at 0(a s ) for m = 2,3,4. However, the resulting expressions are very far from 
compact and at this order we simply choose to use the the information on the moments which 
is available to give us a good estimate of the coefficient function at high x. This confirms that 
again the coefficient function is very peaked for x — > 1 - its size largely compensating for the 
extra power of ais/i^ir). A more sophisticated parameterization than that used in (^J) should 
really include higher powers in ln(l — x), but since the expression matches a range of moments 
very well it will give an accurate representation of the coefficient function convoluted with the 
smooth parton density. 

The predictions for Fl obtained from the parton distributions of the different global fits are 
shown in Fig. 6. The progressive increase at high x is attributable to the large NS coefficient 
functions for x — > 1. At small x the LO and NLOf] predictions mirror the gluon distribution 
(sampled in the region of 2x due to the convolution). The NNLO prediction of Fl also mirrors 
the shape of the gluon at low Q 2 and moderate x, turning over at x ~ 0.05. Then, at even 
smaller x, the very large 0(a s ) contribution of Cl, 9 takes over, which after convolution with 
the gluon, prevents Fl becoming negative and, in fact, results in a steep rise with decreasing 
x. As we evolve up in Q 2 the effect of the 0(a s ) term in Cl,q diminishes and eventually the 
NNLO prediction for Fl mirrors the shape and size of the gluon via the 0(as) term in Cx (ff . 
Hence there is a transition at Q 2 ~ 5 GeV 2 where the NLO overtakes the NNLO prediction^ of 
Fl- At the lowest values of Q 2 the NNLO prediction of Fl should be regarded with caution. If 
we go below Q 2 = 2 GeV 2 the dip in Fl in Fig. 6 becomes negative, indicating the unreliability 
of the NNLO analysis in this domain. 

4 As for P qg , at leading order we present the coefficient of the moment space coefficient function as N — > 0. 
5 Note the very small coefficient of a 2 s /x in Cl, 3 of (10). 



6 ln fact at very low Q 2 and x ~ 10 4 the rate of evolution, dFi/dlnQ 2 , is negative at NNLO. 
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In the region Q > 20 GeV 2 , Fig. 6 shows a LO — > NLO — > NNLO hierarchy in the small 
:c behaviour of Fl, which reflects that observed for the gluon in Fig. 4. As compared to the 
gluon, we see that the NNLO effects in the Fl coefficient function have improved the stability 
of the predictions somewhat. The degree of stability is displayed in Fig. 7, which shows the 
NLO/LO and NNLO/NLO ratios of the F L predictions for two values of Q 2 . The convergence 
is slower for small x, which most likely is due to the influence of missing log(l/a;) terms at 
higher orders. The convergence improves rather slowly with increasing Q 2 . 



6 Predictions for W and Z hadroproduction 

The cross section predictions for W and Z production at the Tevatron and the LHC are shown 
in Fig. 8, together with data from the CDF [|19| and DO f2(| collaborations. The predictions 
labelled LO, NLO and NNLO are defined (schematically) as follows^] 

&lo = ho ® ho 
onlo = hho ® hho ® 1 + as, nlo 

onnlo = /nnlo ® /nnlo ® 1 + «s,nnlo + («5,nnlo) 2 K 1 - 2 ^ (12) 

where the label on«s indicates the order to which the (3— function is evaluated. The NLO and 
NNLO contributions K^ l,2 > are taken from PTJ. The range of NNLO predictions, corresponding 
to the A or B choice for the approximate NNLO splitting functions, is indicated by the width 
of the band. As for Fl, the extrema are given by the AA and BB predictions (see Eq. (2)) with 
the 'average' NNLO partons giving cross sections very close to the centre of the band. Also 
shown in Fig. 8 (as dashed lines) is the 'quasi-NLO' prediction 

Onlo' = /nlo ® /nlo ® 1 + as,NLO K W + («5,nlo) 2 K {2) (13) 



which is the expression used in previous MRST estimates of the W and Z cross sections || |T0|| . 
The NLO' predictions enable us to identify the separate NNLO contributions to the cross 
sections from changing from NLO to NNLO partons and from including the explicit 0(a 2 s ) 
NNLO coefficient functions (K^) in the W, Z cross section perturbation series. 

The LO — > NLO — > NNLO convergence of the predictions is much better than for Fl, because 
the boson cross sections depend mainly on the quark distributions at x ~ 0.05 (Tevatron) and 
x ~ 0.006 (LHC). Since the global fits include high precision F 2 data, there is considerable 
stability in the quark distributions in the sampled x regions, see Fig. 3. 

The jump from cr LO to cr NL o is mainly due to the well-known large 0(ats) double logarithmic 
Drell-Yan K-factor correction arising from soft-gluon emission. The NLO and NNLO cross 
sections are much closer. By comparing with the NLO' predictions, we see that at the Tevatron 



7 A11 quantities are evaluated in the MS factorization and renormalization schemes, with scale choice Q = M- 



v ■ 
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energy the increase of about +4% from NLO to NNLO is due in roughly equal parts to the 
slight increase in the u and d partons in this x range (see Fig. 3), and the net effect of the 
various contributions. 

At the LHC energy the NLO and NNLO predictions are even closer, because (a) the 
contribution is smaller due to an almost complete cancellation between the positive qq and 
negative qg contributions [21|j, and (b) the quark ratios average to unity at x ~ 0.006 for 
Q 2 ~ 10 4 GeV 2 , see Fig. 3. The NNLO band is larger than at the Tevatron because the 
partons are probed at smaller x, where there is more uncertainty in the NNLO evolution. 

We may conclude from Fig. 8 that perturbative convergence is not a dominant uncertainty 
in predicting the W and Z cross sections. This stability indicates the potential value of these 
processes acting as a luminosity monitor for the Tevatron and the LHC. 



7 Conclusions 

In this paper we have taken a first look at a NNLO global parton analysis of deep inelastic 
and related hard scattering data. Although the NNLO splitting functions are not fully known, 
enough information is available to bound their possible behaviour. Even allowing for the full 
spread of the uncertainties of the functions, we are able to draw interesting conclusions. The 
inclusion of NNLO effects gives an overall improvement in the description of the data, which is 
due to the increased scaling violations at both large and small i. In a similar manner, if higher- 
twist contributions are allowed, they decrease in magnitude for both large and small x as we 
increase the order, approaching very small values for x < 0.3, but remaining large and positive 
at large x. The latter behaviour largely reflects the expectations arising from the presence of 
heavy target corrections. 

Fitting to the data using LO, NLO and NNLO frameworks leads to a hierarchy of gluon 
distributions at small x, such that the NNLO (MS) input gluon is found to go negative for 
x < 10~ 3 . However, we stressed that perturbative convergence should be tested for physical 
observables, rather than for the parton distributions themselves. To this end, the LO, NLO 
and NNLO predictions were made for the longitudinal structure function Fl, and for W and 
Z hadroproduction cross-sections. Although the input gluon goes negative for x < 10 -3 , we 
found that Fl is positive for Q 2 > 1 GeV 2 . Despite this the form of the predictions for Fl 
show that the DGLAP approach is not convergent until Q 2 ~ 5 GeV 2 . The convergence 
then improves slowly with increasing Q 2 and reveals a LO — > NLO — > NNLO hierarchy in the 
predictions for Fl, which mirrors that of the gluon but with increased stability. A measure 
of the uncertainty is the ~ 15% change in Fl in going from the NLO to NNLO prediction at 
x ~ 10~ 3 and Q 2 ~ 100 GeV 2 . The convergence deteriorates with decreasing x and most likely 
is due to the neglect of log(l/x) contributions beyond the NNLO DGLAP framework. At low 
Q 2 (Q 2 < 5 GeV 2 ) the log(l/x) terms are even more important. There is also the possibility 
of higher-twist contributions, which for Fl may be different at small x from those for F 2 |27] . 
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On the other hand the predictions of the W and Z hadroproduction cross sections are rather 
stable, due to the more direct relation between the fitted data and the predictions. 

Here we have addressed, in an exploratory fashion, theoretical issues arising from includ- 
ing NNLO corrections in global parton analyses of deep inelastic and related data. However 
new HERA data with increased precision will soon be available. These will be included in a 
new global analysis to yield both an updated set of NLO partons and a first set of NNLO 
distributions. 
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Appendix : NNLO treatment of heavy flavour partons 

For the treatment of heavy flavours we use an approximate NNLO generalization of the Thorne- 



Roberts variable flavour number scheme (VFNS). This scheme was presented in detail in [22 
and the general framework outlined for all orders in perturbation theory. Essentially one obtains 
the VFNS coefficient functions in terms of the fixed flavour number scheme (FFNS) coefficient 
functions and partonic matrix elements A a b. The former are the coefficient functions calculated 
assuming that the heavy quark (denoted by H) has no parton distribution, but may only be 
created via a hard scattering process. The matrix elements define the (jif + l)-flavour parton 
distributions in terms of the n^-flavour parton distributions, i.e. the An a tell one how the 
heavy quark distribution is constructed from the light partons and the A a b,H tell one how the 
light parton distributions are altered by internal heavy quarks (in particular A^} H = 5 a b). The 
VFNS coefficient functions are determined by solving Eqs. (3.5)-(3.9) in the latter of J22|. For 
example, 

n 

r FF(n) _ \ - r VF(n-m) Am) rY VF,PS(n-m) Am) r^VF,NS(n-m) ^VF.PSfn-m}, Am) 

C Hg - 1^ ^Hg ® A gg,H + n ! L 'Hq ® A qg,H + [ C HH + ^HH \ ® A Hg ■ 

m=0 

(14) 

The matrix elements and FFNS coefficient functions are unambiguously calculable, but there 
is some element of choice in the VFNS coefficient functions since there are more degrees of 
freedom than there are constraining equations. One may eliminate this ambiguity by simply 
calculating diagrams assuming one has initial state heavy partons and keeping mass dependent 
terms. However, this leads to unphysical threshold behaviour for the coefficient functions, and 
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we choose instead to impose as physical a constraint as possible. Hence, we make the derivative 
of F^(x, Q 2 ) continuous in the gluon sector (which overwhelmingly dominates) as one switches 
from FFNS to VFNS coefficient functions and turns on the heavy quark parton distribution at 
Q 2 = m 2 H . This choice of VFNS coefficient functions is essentially a freedom in factorization 
schemes, with all schemes becoming identical when summed to all orders, but differing by terms 
~ m 2 H /Q 2 at finite order. 

At NNLO, all VFNSs experience two related technical complications due to internal quark 
loops which may or may not be cut. First, it has long been known that the parton distributions 
become discontinuous at /i 2 = m 2 H at 0(a 2 s ) For example, the heavy quark distribution 

at fi 2 = m 2 H becomes 



f as(m 



(H + H)^ m 2 H ) = {^^)"[A% ® g(m 2 H ) + A% ® q(m%)]. (15) 

The gluon and light quarks also acquire discontinuities as the heavy parton distribution is turned 
on, such that momentum is conserved, see [23|| . These lead to a corresponding discontinuity 
in the coefficient functions, maintaining the continuity of the structure functions, e.g. solving 
flrj) at NNLO at /i 2 = ra\ one obtains 

c £ (2) = c w W+c w S( .,^g, (16) 

The second complication at NNLO arises because the heavy quarks in the final states are no 
longer just those coupling directly to the external vector boson probe, but can be generated 
even when it is a light quark coupling to this probe. In principle it is a technical shortcoming 
of our scheme that the implicit definition of the heavy quark structure function involves the 
heavy quark coupling to the external vector boson. This simplifies the factorization, but is 
not strictly physically correct. A more general prescription is discussed in |[24j| , where a cut 
in invariant mass has to be implemented above which heavy quark-antiquark pairs generated 
away from the external vertex may be defined as observable. 

In this paper we simply ignore both these complications. This is due to the fact that 
the whole analysis is approximate and also because both lead to effects which in practice 
are extremely smalf] - especially when compared to other uncertainties. Both complications 
should be dealt with in a truly precise NNLO analysis once the exact NNLO splitting functions 
are known, though we are confident that they (especially the latter) will lead to tiny effects. 
However, at present we do not even know the NNLO, i.e. 0(a 3 s ), FFNS coefficient functions, 
so a precise VFNS is impossible to define. 

Nevertheless, this is where our heavy flavour prescription comes into its own. Other pre- 
scriptions [F26], |24j which use the coefficient functions from diagrams involving single initial state 



8 This discontinuity begins at 0(as) in some factorization schemes. 

9 The change in parton distributions across threshold was investigated in |^4|, but using GRV98 NLO parton 
distributions (^5|. The discontinuity is dominated by the gluon at small x. The GRV gluon at small x is large 
at n 2 — m 2 , while ours is small and even becomes negative at the same scale. Hence, the effect is very much 
smaller. 
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heavy partons rely on precise cancellations between the heavy quark distributions and terms 
involving the VFNS coefficient functions in order to maintain smooth behaviour. For example, 
at NNLO a large contribution to the heavy quark evolution from a\Pj® needs to be cancelled 
by a term a%P^ \n(^ 2 / m 2 H )C^^} (g) g{^ 2 ) to avoid too quick a growth of F^(x,Q 2 ) for /i 2 
just above m 2 c . In our prescription the correct threshold behaviour is built into C^jjh auto- 
matically, i.e. = if W 2 < 4m 2 ¥ , and such precise cancellations are not necessary - 
simply including NNLO evolution of the heavy quarks without NNLO heavy quark coefficient 
functions at all maintains smooth behaviour. However, we want to obtain the correct NNLO 
high Q 2 limit. Hence, we include the massless 0(a 2 s ) coefficient functions for the heavy quarks, 
but weighted by a factor of f3 — (1 — Am 2 H z/(Q 2 (l — z)) - 5 , i.e. the velocity of the heavy quark 
in the centre of mass system, to impose the correct threshold behaviour at low Q 2 . This pro- 
cedure is very simplistic, but it contains all the relevant physics. Significant improvements to 
this approximate procedure would require the NNLO FFNS coefficient functions. 

Finally, we note that the heavy flavour longitudinal coefficient functions behave like (3 3 , 
and thus are heavily suppressed until very high Q 2 . At such high Q 2 , the 0(a%) coefficient 
functions have become relatively unimportant, and hence we simply omit the 0{a\) longitudinal 
coefficient functions until a more precise analysis is possible. 
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Figure 1: The description of data |12| for the F 2 structure function at a few representative x 
values obtained in the LO, NLO and NNLO global parton analyses. 
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Figure 2: The description of data [12[ for the F 2 structure function at large x obtained in the 
LO, NLO and NNLO global parton analyses. 

17 



Ratio NNLO/NLO 
at Q 2 = 10 GeV 2 




Figure 3: A comparison of partons obtained in the 'central' NNLO analysis with those obtained 
in the NLO fit, first at Q 2 = 10 GeV 2 and then at Q 2 = 10 4 GeV 2 . We show the NNLO/NLO 
ratios for the gluon and the up and down quark distributions. 
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Figure 4: The evolution of the gluon obtained in the LO, NLO and NNLO global analyses. 
The gluons obtained using the extreme forms, A and B, of the NNLO splitting functions are 
shown (dot-dashed curves), together with that from the average (continuous curves). 
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Figure 5: The behaviour of the NNLO contributions to the coefficient function xC^ g {x) for F L 
taking rif — 3. The average of the two extreme behaviours is shown. 
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Figure 6: The predictions for the longitudinal structure function Fl obtained from the LO, 
NLO and NNLO sets of partons. The NNLO extremes (using parton uncertainty alone) are 
plotted only for Q 2 = 2 GeV 2 . 
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Figure 7: The NLO/LO and NNLO/NLO ratios of the predictions of F L , at two different values 
of Q 2 , shown to indicate the degree of perturbative stability of the analysis. 
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Figure 8: The predictions of the cross sections for W and Z production and leptonic decay 
at the Tevatron and the LHC obtained from parton sets of the LO, NLO and NNLO global 
analyses. The cross sections labelled LO, NLO, NLO' (dashed line) and NNLO are defined in 
Eqs. ( ^2|JT3| ). The band of NNLO predictions corresponds to the A, B variation of the small- 
er approximate splitting functions, as discussed in the text. Also shown are measurements 
obtained at the Tevatron 20fl . We take the leptonic branching ratios B{W — > Iv) = 0.1084 
and B(Z -> l+l~) = 0.03364. 

23 



